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INTRODUCTION 



The irregular motions, called turbulence, play a 
prominent part' in all technically important ^low phenom- 
ena. Turbulence, on the one hand, is the cause of unde- 
sirable flow resistance, while, on the other hand, it has 
the very useful.' characteristic of increasing the pressure 
in the currents. The control of these phenomena is very 
important for the flow specialist. Nutnerous researches 
have therefore been recently undertaken for the purpose 
of discovering the law's; of turbulent flow. In the present 
article an attempt is' made to' review the most important re- 
sults of these researches. Relations of immediate practi- 
cal interest are discussed. 

The first two^ sections treat of two prominent ques- 
tions, namely the origin of turbulence and the character- 
istics of turbulent currents. In the third section con- 
clusions are drawn for the flow along a rough wall, where- 
by an important relation for the veiocity distribution is 
revealed. The principles are also applied to- straight 
rough and smooth tubes. Here it was possible to develop 
formulas for flow velocity anti resistance, which show ex- 
cellent agreement with the experiments, and which also in 
contrast with previous purely empirical formulas, hold 
good for very large Reynolds Numbers for which no experi- 
mental data are available. The peculiarities in tubes 
with fine-grained roughness at moderate Reynolds Numbers 
are represente'd by a single curve. Test results with arti- 
ficially roughened tubes are given and confirm the relation- 
ship mentioned. 



^'Neuere Ergebnisse der Turbul enzf or schung . " Zeitschrift 
Qes Vereinos deutscher Ingenieure, vol.. 7, no. 5, February 
4, 1933, pp. 105-114. 



2 IT.A.'C.A. Technical Memorandum No. 720 



The resnlts ottained with tubes are applied to the 
resistance of plates to a longitudinal flow. Moreover, 
the characteristics of the flow in wide and narrow and 
curved channels, as likewise the mixture phenomena of 
fluid currents with surrounding fluids and also the phe- 
nomena behind moving bodies are considered. Lastly, newly 
discovered relations between the turbulent exchange of ve- 
locity and heat are considered, and new conclusions are 
drawn regarding the finer details of turbulent flow. 

Durin^:; the last decade the investigation of the ir- . 
regular mixing motions, which are called turbulence and 
which affect all technically important flows, have been 
especially thorough and fruitful • These mixing motions 
produce effects, as if the viscosity of the fluid were in- 
creased a hundred or ten thousand fold or even more. This 
circumstance causes the great resistance of fluids in 
pipes, the frictional resistance of ships and airships and 
other undesirable resistances, but also the possibility of 
increased pressure in diffusers or along airplane wings 
and blower vanes. Without turbulence, separation wpuld oc- 
cur in these cases, so that there would be only a small re- 
covery of energy in the diffuser and impaired efficiency 
of wings or vanes. 

The investigation consisted of a determination of the 
numerical data and their systematic arrangement. General- 
ly the investigation was not carried to an actual theory 
(which is very difficult), but the results help to support 
theoretical conclusions. Often dimensionsal considerations 
together with intuitive insight lead to important conclu- 
sions. If, e.g. density (i.o. inertia) and viscosity are 
the only determinative properties of the fluid for the phe- 
nomenon, one is led to a Eeynods Number = density/viscosity 
X velocity X length (Ee = vl/y, in which is the "l-ine- 

matic viscosity", i.e. viscosity/density). If Reynolds 
Number has the same numerical value in two cases, we may 
expect exactly the same course in both cases, only with a 
different length and time scale according t o c ircums tanc es . 
In individual cases the application of this rule may, how- 
ever, require consideration as to Which velocity and 
which length is actually determinative for the process. 

There are tTO main questions which were investigated 
theoretically and experimentally : 

1. How and under what immediate conditions does tur- 
bulence originate? 
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2. What can he said regardinf^ t\irh\Tlent motion, par>- 
ticularly regarding the ^aean valne.g of the velocities and 
f or ce s ? 

. The second question is ohvi.ously t;ie more important 
one from the technical viewpoint. 



ORIGIN OF TURBULENCE 

• - ■ ' ■ . . • ■ - . 

^ Regarding the first question I can be ,quite hrief, 

both because I have recently expressed myself on this sub- 
ject (reference 5) and because there is here much that is 
still in doubt. The most important fact is that turbulence 
always occurs when the velocity profile shows, a turning 
point (fig. 1) and when the viscosity effects are not too 

^ great. Any flow with such a velocity profile is unstable 

in the absence of fluid friction, i. e. small deviations 
in magnitude and^ Air ec t ion increase of themselves and cause 
a complete reversal of the flow. An originally slight wave 
in the streamlines leads .gradually to the production of 
turbulence .through tlie toppling ov,.er of . the . waves . These 
p.henomena can be delayed, by. strong viscosity effects. 

This indicates that the tendency to become tubulent 
depends on the .magnitude of. the Reynolds- Numb.er. Velocity 
profiles with turning point .occur, e.g. in the bou-ndary 
layers produced by vis.cosity effects, when th-e. pressure in- 
creases in. the directio.n of flow or, in other words, when 

^ the flow is retarded. Such points in the fluid therefore. 

have a strong tendency to become turbulent, but even the 
unaccelerated rectilinear flow along a wall tends to become 
tu^rbulent- at a sufficiently large Reynolds Number. This 
can bo explained by the fact that the inflow is never abso- 
lutely undisturbed and that there are always, some irregu- 
larities in the velocity distribution. Unstable velocity 
distribution is largely' due to only slightly damped turn- 
■ ing motions with axes parallel to the direction of flow. 
Such turning motions direct some portions of the fluid a- 
gainst the wall and other portions away from it, so that, 
even at low velocity, with the lapse of time, pcrtipns hav- 

^ ing a lower velocity become interspersed with portions hav- 

ing a higher velocity, thus rifecessar ily produc ing instabil- 
ity. 

' There is still another cause of turbulence, which was 

* 4 



4 
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discovered in a theoretical manner (references 3, 3, and 
4) and which call for special consideration when there 
are none of the ahove-ment ioned dis tur.tanc es . In the flow 
along a wall, there occur certain slow d is t rii"banc es which 
aoove a certain critical Heynolds K"amher, increase in 
strength and thus produce in their retarded zones, after 
their amplitudes have "become great enough, the preliminary 
condition for turbulence. It is worthy of note that the 
critical Reynolds Numbers- for two different cases, as de- 
termined theoretically by Tollnien (reference 4) and 
Schlichting ( reference 6) , are in good agreement with the 
experimental values. 

Experiments on the Pr oduc t ion of Turbulence 

In order to obtain more light on this question, we 
investigated the production of turbulence by experiments 
in channels 20 cm (7.87 in.) wide and 6 m (19.68 ft.) long. 
Though we proceeded with great dare, we found it impossible 
to eliminate all the disturbances, so that here and there 
nuclei of turbulent motion developed in irregular succes- 
sion and spread quite rapidly. 

Clearer pictures were obtained by purposely initiating 
a disturbance in the flow, as, e.g. by adding or removing 
a little water through a small piece of screening inserted 
in the wall. In the first case, when a small amount of 
water, not yet participating in the flow, is thrust between 
the wall and the moving mass, instability is immediately 
produced and turbulence develops at the point of entrance. 
The amount of water introduced may be very small. In the 
second case the greatest disturbance occurred in the por- 
tion of the flowing water opposite the screen at the be- 
ginning of the removal by suction. Behind this point the 
thickness of the boundary layer was reduced by the suction, 
and the inner portion of the water flowing past the bound- 
ary layer had therefore to flow over a sort of step from 
the thinner boundary layer to the thicker layer on the 
'dov/nstream side of the screen. This created enough of a 
disturbance to cause the disintegration of the boundary 
layer in a short time. Figure 2 shows this effect and the 
further develop?'ient of the turbulent region. 

The flow was rendered visible by scattering aluminum 
dust on the surface of the water. A slowly operating mo- 
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t ion-picture camera was monnted on a car which kept pace 
with the flow so that the same group of vortices remained 
in the field of the camera. In the top picture the oh- 
lique streamlines at the left show the location of the suc- 
tion point, while the formation of the first vortex in the 
middle indicates the location of the "step./ Other vorti- 
ces developed on the upstream side. In the last picture 
the original vortex is shown at the extreme right. It is 
evident that - it carried water from the "boundary layer 
(which was purposely strewn more thickly with aluminum 
dust) far into the interior of the flow. 

CHARACTERISTICS OF TURBULEJfT CURRE5TTS 

We will now consider the l^ws of fully developed tut- 
hulence. The method of presentation which I shall employ, 
does not follow the historical, development, btxt is intend- 
ed to show the present status all the more palinly. I 
shall hegin with a statement, regarding the behavior of an 
ideal fluid without viscosity. In reality there is no 
such fluid, but it is of advantage for many considerations 
to know what would occur in such an ideal fluid, because 
the laws of the .ideal fluid (due to .the absence of viscos- 
ity) are sim.pler than those of an actual fluid. 

According to our previous statements the tendency to 
creases or, in other words, as the viscosity decreases 
(under otherwise like conditions). At the zero limit of 
viscosity -the Reynolds Number obviously becomes infinite, 
necessitating the c one lus ion that the flow of an ideal . 
fluid, would generally be turbulent. If it is also, assumed 
that the bodies or walls, past which the fluid flows, are 
mathematically smooth, the surface friction would also be 
zero and we would thus obtain the theoretical behavior of 
the ideal fluid, as stated in old textbooks on hydrodynam- 
ics. If however, the siirfaces are rough, it m.ay be assumed 
that an area of separation develops at each individual 
point of roughness, however slight.*'^ The flow thus acquire 

^In slightly viscous fluids, a regular separation of the 
flow occurs on projecting parts o'f' the. wall. In the lim.it- 
ing transition to vanishin^^; viscosity,. Helmholtz separation 
surfaces a re developed with^finito velocity jumps. 
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a turbulent charac t er .f rom the mutual effect of the various 
sm"all areas of separation which are unstable in themselves 
and have a disturbing effect on one another* , At. each rough 
spot a pressure difference develops between. its upstream 
and downstream sides, thus producing a resistance which is 
proportional to the square of the velocity. 

From this consideration it may be assumed that it is 
permissible to make theoretical assumptions regarding the 
laws of turbulence, in which the viscosity of the fluid is 
put at zero. The following considerations clearly show 
that we are thus on the right track and that, as a matter 
of fact, the turbulent resistance in the interior of the 
flow is practically independent of the viscosity. In a 
thin layer near the wall, however, the effect of the vis- 
cosity persists, provided it- is not concealed by the ef- 
fect of great roughness. 

We will briefly explain a conception which has been 
found'useful for the more accurate investigation of the 
turbulent mixing processes. This is the so-called "mix- 
ing path," which plays a similar role in turbulent mixing 
processes to that played by the mean free path in the mo- 
lecular diffiasion of gases. In both these precesses shear- 
ing stresses (or apparent shearing stresses) are developed 
by the continuous interchange of energy between fluid lay- 
ers flowing parallel to one another at different velocities. 
The following simplified representation can be made of 
these really quite complex processes. 

It is assumed that any particle, which, by collision 
with neighboring particles, acquires a motion crosswise to 
the flow, has, in the direction of flow, the mean momentum 
of the layer from which it came, and that it now traverses 
a distance I crosswise to the flow, before it collides 
with other particles or mingles with them. Such exchanges 
occur in both directions, and thus the faster layer re- 
ceives particles from the slower layer, which naturally re- 
tard the former, and, conversely , the slower layer receives 
particles from the faster layer with an accelerating effect 
on the former. 

The effect of the two fluid Layers on each ot.hBr is 
therefore the same as if there were friction between them. 
The difference between the nolec":;lar processes and the 
turbulent processes is due only to the fact that, in one 
case, the individual moleciiies, and, in the other case, 
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whole groMjjS of molecules participate in the exchange. If 
u is the velocity of the flow and y the coordinate in 
the direction at right angles to the flow in which the 
chsinge in velocity occurs, the difference hetween the veloc- 
ities of the • two layers, separated by the distance I, is 
I (du/dy) . This, according to what precedes, is also the 
velocity difference of a particle which, coming, from the 
other layer, mingles anew with its present environment, 

' In order to determine the magnitude of the frictional 
force or, more accurately stated, the shearing stress "be- 
tween the two layers, we must know the magnitude of the 
mass exchanged per second. This, as referred to the unit 
area, can be expressed by the product of the density 
pC~7/"s) and an exchange velocity v'. In the case of the 
molecular motion, this velocity is proportional to the ve- 
locity of heat transfer. Since the latter, is one third 
each along the x, y and z axes and since, in our oxampl^e, 
re can put, in first approximation*, v' - c/3, where c 
is the mean velocity of the heat transfer. Hence the shear- 
ing stress*^ 

1 , du ^ du • /t \ 

T = -r p c I T- = T^, J- (1 ) 

dy 'dy - 

In the case of the turbulent exchange of masses, the 
velocity v* should naturally be taken of the same order 
of magnitude as the difference in t he veloc i t ies of the 
two layers at the distance I from each other, since the 
fluid masses collide at volocitic^s of this order of mag- 
nitude (references 8, 9, and 10). On eliminating the un- 
known numerical factor v», we thus obtain the shearing 
stress • 



The elimination of the numerical factor only denotes a 
somewhat different definition of I . In this way wo ob- 
tain, for the simple viscosity effect, shearing stresses 
proportional to du/dy- and, for the turbulent exchange 
(whereby the effect of viscosity is disregarded), shearing 
stresses proportional to (du/dy) ^, which is in good a- 
gr e enen t _^i 5.^_t h e _hy d raul _r e s i stances proportional to the 

By a more accurate calculatiqn, Boltzmann found, for ther- 
viscosity ri , the value r] = 0.3503 pel. which differs 
but little from that inequation 1. 
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square of the velocity. 

With formula 2 the problem of the hydraulic flow re- 
sistances is "brought "bad: to the other problem of the dis- 
tribution of the mixing path I in the flow. So long as 
we have no rational theory of turbulent flow which deduces 
the laws of turbulent phenomena from hydrodynamic differ- 
ential equations, we have to obtain the data regarding the 
distribution of the mixing path by experimentation, so that 
only one unknown quantity is thus replaced by another, Nev- 
ertheless, considerable progress has been made, since it 
has been found, at least for the larger Reynolds Numbers 
(from about 10^ up) that the mixing path is practically 
independent of the magnitude of the velocity and is, more- 
over, subject to quite simple rules for its distribution 
in space. 

Dimensional considerations often furnish useful in- 
dications. For example, in considering the flow near a 
0 r / 

more/less smooth flat wall, on the assumption th^.t neither 
the viscosity nor the roughness/of the wall has any appre- 
ciable effect at the point under consideration in the in- 
terior of the fluid, we are in a position to make a state- 
ment regarding the distribution along the mixing path. For 
a point at the distance y from the wall there is no other 
characteristic length than this distance y. The mixing 
path I is also a length, so that there is no other pos- 
sibility than to put the mixing path proportional to the 
distance from the wall: 

I ■ = Ky. 

Here k is a universal numerical 
be determined experimentally. If 
flow in which the shearing stress 
tain 



d u ^ ^ I 
dy *^y a/ p 

according to equation 2, and therefore 

u = ^ 7 ~ (In y .+ const . ) (3) 

Such a velocity curve, depev^dent on the distance from the 
wall, is quite li'ze the one actually observed (fig. 3). 



coefficient, which can 
we assume a state of 
T is constant, we ob- 
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Comparison with the experimental resrats yields the n^jmher 
0.4 as the approximate va,lv.e of k. 



barman's Theory 

Von Karman (reference 12) assumed that the turbulent 
mixing processes are the same in all cases, so that only 
variations in the len:5th and tim.v scales, occur from case 
to case and from place to place in the flow. Under these 
circumstances the effects of viscosity are regarded as neg- 
ligible in cornpariaor with the effects of turbulence. Con- 
clusions are now drawn from Euler^s equations regarding 
these two scales, the first of which obviously agrees in" 
principle with our mixing path I. The velocity u of the 
basic flow, which is assumed to be a function of y alone, 
is determined from a Taylor series interrupted after the 
quadratic term. The :r:ean forward velocity of the x^rticle 
under consideration has no immediate effect on its inner 
motion. Of the given quantities therefore, only du/dy and 
d-'-u/dy2 need to be considered here. We first have a time 

• T ~ 1/ 

oy 

as the time criterion for tho period of the mixing process. 
For dimensional reasons, the interference velocities u* 
in the X direction and v< . in the Y direction are therefore 
proportional to I/t, i.e. 

da 
-oy 

which agrees with the previous formulas. For the longitu- 
dinal scale of the mixing process, Von Karrnan finds thu re- 
lation 

dy' dy*^ 

in which is a constant det ermined experimentally . 

This expression of Karnan's theory go^s beyond previous 
expressions, because it furnishes a method for calculating 
the magnit^ide of the mixing patn independently of the dis- 
tance from a wall. If this expression is :lnt rod^ic ed into 
equation 2 and integrated on the assumption of a constant 
shearing str^vss in the rf:,gion under consideration, we ob- 
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tain 

i.e- practically equation 3 again. The required agreement 
with the experimental results ohviously leads to putting 

K. Hence both formulas yield the same velocity dis- 
tribution in case of constant shearing stress. 

There is no longer any agreement regarding the shear- 
ing stress in the other assumptions. lioreover, the formu- 
la I = Ky is yrithout any valid basis, since, due to the 
variability of the shearing stress, a still further length 

jjf^l is available; but even Karnan's formula 

du/dfu 
^ - ^- dy^d y''^ 

here means only another estimated approximation,^ since it 
was obtained by di sr e^^ir ding the sffect of d^u/dy^ and 
higher terms in the ^ezres 'iov el :;pment for u. In the 
case T = cOii^tart, the twc solution*^, coincide, because 
the velocity distribution, according to equation 3, is trans 
f erred by changing the inte^'ration constant, in case the 
shearing str3>:«s t r'c^mains unaltered, so that there is 
here also a pronounced similarity with the basic flow. 

From equation 3 it is easily seen that the quantity 
^t/p is a velocity. This velocity is very valuable for^ 
various similarity considerations in what follows. We will 
therefore designate it by v^ and call it ^' shear ing- st ress 
velocity." The formula T = p v* ^ is of similar form to 
that far the dynamic pressure 

which is comprehensible for dimensional reasons, sinre the 
shearing stress is also a force per unit area. The appar- 
ent shearing stress T of the turbulence is generally very 
small as compared with the dynamic pressure. Hence in 
we are also dealing with a velocity \vhich is relatively 
small as compared with the flow velocity u. Comparison 
with equation 2 sho\7S, moreover, that 
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Hence v,^ is of the order of magnitiide of the mixing ve- 
1-Ccities. ••■a* - an;d • v • ; - 



^^LOW ALOHG A ROUGH WALL 



From our standpoint the flow along a rough wall is 
simpler than along a smooth wall, "because the viscosity 
plays a preponderant role in the latter case, hut not in' 
the former. It is theref ore bet t er to consider the flow 
along a rough wall first. If.', k is a length indicating 
the roughness of '.the wall, it", follows, from a simple sim- 
ilarity consideration on the^'hasis of the id^al fluid,, 
that the velocity distributions near the wall, with geo- 
metrically similar roughnesses, are also geometrically sim- 
ilar, so that the si7.e of the grain k furnishes the cri- 
terion for it. The formulated expression of this relation- 
ship is that the velocity at the distance y is a function 
of the ratio y/k. If this velocity distribution is hased 
on equation 3, which, according to what has preceded, is at 
least advisable for the regions farther in the interior of 
the fluid, it is found that the integration constant of e- 
quation 3 = constant - In k. 

A hitherto unpublished series of experiments by Niku- 
radse with tubes of various diameters, which were giv-3n 
different degrees of roughness by gluing to them sifted 
sand with a suitable varnish, showed that the new constant 
- 3.4 In 30, k being the mean diameter of the grains 
of sand used to oroducc the roughness. With l/^= 2.5, 
we obtain the formula 

u = a.5 v^lnf^r-^^ 

By a shifting of the coordinates by the amount of k/30, 
it is also possible to obtain u = 0 for y = 0.'^ rTonce, 

30 V \ 

u = 3.5 v^ln yl ~— y ^ ^'^ 

"^It is still an open, thou.f:;h not. very important question as 
to the exact location of the axis of t;he coordinates between 
the protuberances of the roughness. 



12 
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or, if the natural logarithm is replaced hy a common log- 
arithm, 

u = 5.75 log (1 + — ) (5a) 

Equations 5 and 5a therefore show a fixed relation between 
the velocity distribution, shear ing- st r e s s velocity, dis- 
tance from the 77all and the degree of roughness k. This 
first holds good for the kinds of roughness used in the 
experiments. For other forms of surface roughness, more- 
over, there is prohahly another . number instead of 30, also 
dependent on the manner of defining the roughness scale. 
Preparations for tests in this connection are being made 
in Gbttingen. 

Equation 5 immediately affords us the opportunity to 
check the above statement regarding the behavior of the 
ideal fluid. Represent the velocity at the distance y = h 
by u u^. With this assumption v,^ can be eliminated 
from equation 5a. 



and consequently 



u 



5.75 log 



1 + 30 > 



k 



log^l + 30 



u ~ u 

' log(l + 30 



(6) 



The corresponding shearing stress is 

P ^i' 



T = p v^ = 

S3 



log (^1 + 30 |j] 



(7) 



from which it follows that the shearing stress is propor- 
tional to the square of the flow velocity u . The effect 
of the roughness of the wall is likewise shown by equation 
7 . 

If we pass to the mat;iema t ical ly smooth wall, i.e., 
to k = 0, then, according to equation 6, u = u^^ and 
T = 0 for all valTies of y constant, as stated in the 
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classical hydrodynamics on the ideal fluid. It is also 
obvious that even a submicroscopic roughness with a k 
of the order of ma{^;nitude of the diameter of an atom would 
still show considerable deviations from the ideal behavior. 
0\iv formulas can.no longer be used for such cases. The 
relations are here considerably altered by the viscosity, 
as will be shown in what follows. 



THE FLOW IN TUBES 



It is an important discovery that, in a straight tube, 
the relative motion of the fluid particles at moderately 
large Reynolds Numbers depends on the fall in pressure and 
not at all on the character of the wall, so that therefore, 
with constant fall in pressure, the velocity-distribution 
curves in tubes of greater and less wall roughness can be 
brought into conformity by shifting along the velocity ax- 
is (of course asic?e from a layer in immediate contact with 
the wall, where the velocity increase is naturally greater 
on a smoother surface than on a rougher one). This rela- 
tion was discovered byDarcy (reference 14) 75 years ago 
in his researches . on resistance in pipes and was then em- 
phasized, but was subsequently forgotten. Fritsch dis- 
covered it anew by direct observation in his experiments 
with rough channels at Aachen {reference 17). Prom our 
standpoint this discovery is identical with the fact that 
the distribution of the mixing path along the inside of the 
tube is practically independent of the nature of the wall. 
In connection with our earlier discoveries it is natural to 
S'.irmise that the formula 



I = r f^^ 



can be written for the mixing path, where y is the dis- 
tance from the wall and r the radius of the tube. Since 
the distribution of the shearing stress along the tube is 
known when the pressure fall is given, the distribution of 
I can be verified by measuring the velocity distribution 
with the aid of eq-'^ation 2. It is found that the above 
stateinent is confirmed, at least for the hi/^her Reynolds 
Numbers. Fiijpj.re 4 gives the result in nondimons ional form 
and consequently shows the course of the function f^ . In 
the function f.^(y/r), l is the m.ixing path; r, radius 
of tube; y distance from wall; and k, miean longitudinal 
dimension of toughness. 



14 
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Conversely, on the "basis of this function and with 
the aid of equation 2, we can calculate dn/dy, from 
which, "by an integration, an expression for the velocity 
itself can he obtained. On the introduction of the shear- 
ing stress velocity v* this expression takes the form 



V, 



max 



* 2 vr / 



(8) 



This equation, which was first developed hy Von Karman 
(reference 12), has also "been experimentally confirmed, 
as shown hy figure 5, in which the test points are given 
for smooth tuhes and for various rough tuT-es. In function 
fgCy/r), "^max ^-^^^ maximum flow velocity; u, flow ve- 
locity at the point y; v* , shearing- st res s distrioution • 

shearing stress; p, , density. 



We can now pass from the velocity u at any distance 
y from the wall to the mean velocity u. We thus ohtain 
from equation 8 an expression of the form 



u - u 

max 



X coefficient 



(9) 



Nikuradse's Gottingen experiments yielded 4.07 as the val- 
ue of this coefficient. It was a piece of good luck that 
our equation 3 or the special form for a rough wall (equa- 
tion 5) yielded^, up to the middle of the tube, a useful 
approximation for the function t^{y/r), namely 

5For more accurate calculations, a small supplementary terr>.i 
would have to he added, which will he included later, at 
least in the final result. 

^Darcy (reference 14) deduces from his experiments 



^max - ^ = 11-3 



(r 



y) 



3/: 



( 



is the gradient and therefore = - gp dx' meter is 



the unit of length) • This equation can he put in the form 
of equation 8 and thus becomes 



f ^ = b.OB {' 1 - ^ ^ 
2 r V J 



which, with the exception of the wall vicinity where Darcy 
made no measuirement s , a;":roes .very well with modern results 
(fig. 5). 
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^ f i^'^, - 2.5 m - 5.75 log J (10) 

We npw have all that is needed to calculate the resistance 
of a tcagh tn"be for a given quantity. We will first write 
the customary expression for the drag coefficient \: 

^ dp _ \ PU^' / n . N 

dx, 6. "2 

F.rbm the e qu il ihr iuni of a water cylinder of radius r = d/2 
we ohtkin, for th'e shearing stress t • of" the wall> the 
Express ion 

' ■ ^" " - dp . ^ 

- TT -r~- ^ 2 TT r T • 
ax 

and acc ordin:^rly 

dx r r 

Th'e comparison of equations 11 and 12 yields, with the 
tube diameter d = 2 r, 

■ ^ ^"^^ 

By .the .u§e ,Qf equation 5a at the middle of the tuhe (y ~ r) 
we. obtain, .when, \inder the logarithm, we disregard 1 in 
comparison with the very great value 30 r/k and put log 
30 = . 1 .477 , 

-^max '-^ v^(5.75 log ^ + 8.5) ' (.14) 

On the other hand, according to equation 9 

U = Ui^ax ^ "^'^^ " v^(5.75 log ~ 4- 4.43) 
T;^king equation 13 into consideration, we now have 

\ = = ^ — • (15) 

(5.75 log I -f 4.43)"^ (2.0 log 1+ 1.57)" 

k ■' ^ k ■ 



This is very well confirmed by experiment, with only the 
slight difference that 1.74 is better than' 1.57 in the de- 
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nominator. This difference is connected with the sup- 
pressed auxiliary term in equation 10. The experimental 
conf irmat_ion of the formula is hest accomplished hy plot- 
. t ing ; l/ /v/X against Ipg r/k. According to the foregoing 

' i/a/X = -2.0 log I 1.74 (16) 

The plotting must therefore yield a straight line. Figure 
6 shoYTS this line for six rough tubes according to measure- 
ments "by Nikuradse. (See ^also figure 9.) The general form 
of eqiaation 14, as likewise an equation analogous to equa- ^ 
tion 16 for a coefficient of resistance based on ^j^^j^^* 
first developed by Von Karman. He also made the rectilin- 
ear graph. 



Effect of Viscosity (smooth tube) 



It has already been mentioned that the effect of vis- 
cosity is greater when the roughness is less, but of course 
only on the boundary- lay er phenomena. The rough places are 
here more or less covered by a slower-moving layer of fluid 
and are thus rendered ineffective as regards resistance. 
Progress can also be made here with a dimensional consider- 
ation. The shearing stress is responsible for what takes 
place on the wall and consequently the velocity v* based 
on this shearing stress, and also the criterion of rough- 
ness k. A wall characteristic v^k/i5 can be developed 
from these two with the kinematic viscosity by analogy with 
the Heynolds Number. Since, with fixed v* , the state of 
flow in the interior remains unaltered, the only remaining 
problem is to adapt the integration constant of equation 3 
to the new relations. This is accomplished by introducing 
a modified roughness criterion. 



k» = k f 



/vi> k N 



3V ^ 



instead of k, into equations 5 to 7 and 14 to 16. Regard- 
ing the course of the function fi,, it follows from the 
foregoing that it must be equal to 1 for large values of 
the wall characteristic, in order to restore the previous 
relations. An immediate conclusion can, however, be drawn 



as to what 

values of 



form the function 



must assume for small 



•The observations show that, for slight 
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but still appreciable roiaghness, the rough tube does not 
dif f er pr^ct ica], ly from a perfectly sm.ooth tu^^t provided 
theyEeynoid^ , Ihiniber . 1 s not umisually high^ . 'Such a .CQrldi- 
t ion . "i s ". ob tai.n ed when . /" \ 

■ • • f M ^^-^0 coefficient x^^^-- 

since k is thus removed from the foregoing formulas and 
is replaced by coefficient X The experiments confirm 

• V^ . • : 

thif5 res"'Tlt and show, w i th - r e sp. e.c t to, the coefficient 
which leaves the dimensional consideration still open, 
that our previous value of k/30 must be replaced by 
v/9 . Instead of equation 5a, we now obtain the formula 
for the velocity distribution in the tube 

^ = (5.?5 log 4- 5.5) (1?) 

On plotting u/v^ against log v^y/i" , we obtain a straight 
line which must contain all the points near the wall for 
the velocity profiles of all smooth pipes. An exception is 
formed only by the values at very small nond iraens ional dis- 
tances from the wall y^y/v, at which the turbulence is 
still affected by the viscosity. Up t o the previously 
mentioned supplementary function, equation 17 is also valid 
to the middle of the tube. The experimental points in fig- 
ure 7 actually contain not only the parts near the wall* 
but extend almost to the middle of the tube. One can there- 
fore note small systematic deviations from the straight 
line, which of course have to be considered in a more ac- 
curate theory (reference 18), 

For comparison figure 7 also shows, by a dash line, 
the veloc it j-di str ibut ion law 



. 8:7ll^S' (18) 



A 



'1 / 7 



/ 



as determined on the basis of the Blasius formula for the 
friction of the tube. It is found that, in a central re-- 
gion for which alone data were formerly available, it prac- 
tically coincides with tne straight line of equation 17, 
but deviates considerably abovo and below this region. In 
fact it was long since discovered that, at higher Reynolds 
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Num'bers, the seventh root is replaced "by the eighth and 
nixith roo t s., etc . The' refiL'Son ' f or this behavior is mani- 
fest s ince the law of the seventh root now appears to he 
only an approxima t ion f o rmula for the real law, which is 
represented hy equation 17, wherehy the particular numeri- 
cal values of the approximation formula naturally still 
depend on the region in which they should agree with the 
accurate formula''. 



For the coefficient of resi stance , we obtain from 
equation 16 by the same modification 

1 r 
=2.0 log — +0.5 

Taking equation 13 into consideration, we can put 
v^ r _ u r v^ u d 1 



V 11 V 2 J 8 

With u d/v = Re, we obtain 

l/JY 2.0 log (Re^/TD - 1.0 (19) 



This formula was verified experimentally be Nikuradse 
(reference 20) up to the Reynolds Number 3.4 X 10^. It 
must be changed only by the consideration of the previously 
mentioned supplementary function of the numerical value 
from - 1.0 to - 0.8. The final formula for the resistance 
coefficient is then 



l//K = 2.0 log (Re ^/X) - 0.8 (20) 

The calculation of the resistance coefficient corresponding 
to any given value of Reynolds Number encounters no partic- 
ular difficulties, although ^/ K occurs once more on the 

right side. _0ne can, for example, assume provisionally any - 
value for on the right side and calculate 1 ^/X and 

then repeat the process, if the discrepancy is too great. 
In figure 8 the course of X is plotted with respect to Re 
according to equation 20 together with the experimental val- 

^Sc^iow log v^y/ v= 2, the .straight line of equation 17 
shows appreciable deviations from the test, points. This is 
due to the influence of the viscosity. If the. smallest su- 
percritical Reynolds Numbers are disregarded, this deviation 
occurs only in a very thin layer near the wall of the tube. 
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ues. By especially good iMck this formula agrees with the 
exporimonts doTzn to the srunllest sup'^rcr i t i cal Reynolds 
Numhers. 

We now turn once more to the general prohlem of the 
rongh tuhe. On the "basis of iriea sur erient s by rlikijraclse 
(now "being prepared for pul'l icat ion ) the course of the re- 
sistance coefficient is plotted in figure 9 against the 
Reynolds N\im"ber for tuhes of different relative roughness 
k/r. The curves in figure 9 are "based on experiments with 
tu"bes of well-defined roughness produced "by gluing grains 
of sand of definite and different sizes (k) to the inside 
of tu'bes. The conditions to the left of the critical Rey- 
nolds Nurn"ber represent the laminar condition of smooth 
flow. It is evident' that th.ere is hero very little differ- 
ence "between the smooth and rough tu"bes. The curves di- 
verge greatly, however, as soon as the tur'bulence "begins, 
i.e. a"bove Re^^^^^ The curves for the lesser roiighness 
first follow the curve for the smooth tuhe and then sepa- 
rate from the latter in order. 

The foregoing considerations indicate a way to find 
a law for the tur"Du.lent portion. We will take the wall 
characteristic v^k/i/ or its logarithm as the; f.ibsc is sa and 
a quantity which is constant accordin;;; to the laws of the 
fully developed rou<^rhness flow as the ordinate. For exam- 
ple, we can take the quantity 

- 2.0 log ^ 

or, if .we want the corresponding law for the velocity dis- 
tri"bution» the quantity 

- - 5.75 log 

V 'r 

The plotting of these two quantities on the "basis of fiio 
experimental results "brings in fact the test. points meas- 
ured with very different roughnesses approximately on a 
single curve. T'he two ciirvts agreu with each other up to 
the scale corr esp)onding to tht relations here reprensent ed. 
The whole pro"b'lem thus finds a very comprehensive solution 
on the tpsis of combining a f ew experimental values with 
theoretical conclusions. What remains to be done is to 
find the curves for other forms of rou.^hness in addition to 
the curve of figure 10, which we have thus far determined 
only for the special sandpaper form of roughness. Prepara- 
tions are now being made for sucn experiments. 
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APPLICJ^TION, TO OTHER CASES 
Plate Resistance - Accelerated and Retarded Flows 



From the behavior of the flow in tubes, when the 
P'lasitis law of resistance 



\ = 0 .316 



-1/4 

/u d\ 



dominated the field, concl-as 5. ohs had already "been drawn 
regardin/^ the frictional resist<ance of plates suhjected 
to fltiw along their surface (references 15 and 16). Ac- 
cording to the momentum theory, the decrease in the mo- 
mentum of the flow due to the friction was represented 
"by a formula in terms of the exposed length of the plate 
in accord with the' laws for the velocity distribution. 
This decrease in momentum per unit. len^^:th alont^:^. tho plate 
was expressed as equal to -the frictional force per unit 
len^-th. The resulting formula for the coefficient of 
frictional resistance c^ (resistance divided "by the sur- 
face area and dynamic pressure) ; 

= 0.074 fl-iV^'' (21) 

(I- length of plate, v = velocity of plate), showed simi- 
lar discrepancies, in comparison with the experimental re- 
sults, to those shown in the resistance of tu"bes. The oh- 
vious thing to do now was to apply the improved law of tubu- 
lar flow also to plates. The calculations are here rather 
troublesome. They were first made by Von Karman (references 
13 and 2l) . A new calculation in a somewhat different ^:ay 
ras made by the writer ( reference . 18 ) , who compiled a nu- 
merical table the values of which agree very satisfactorily 
with Kempf's measurements. The values in the table were 
obtained by the . f ollowing approximation formula of H. 
Schlichting, which though it- is only an interpolation for- 
itiula, cen 06 used throughout the whole practical region of 
turbulent, flow. 

= = (22) 

1 0 1: ~ i . ' 



Por.the rongh plate, a corresponding calculation was made 
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OD tire "basis' of the law of rou^'^ihness represented in figure 
IC (reference 19). 

The "behavior of the tiir"b.j.l ent friction layer in an ac- 
celerated or retarded flow is of greater importance. An 
important special case, the flow in a widened or narrowed 
channel with flat side walls, was investigated "by Bench for 
air (reference 22) and "by Nikuradse for water (reference 
23). I^nri's v^orlc at Znrich sho-ild "be metioned here, as al- 
so C-^ino's experiments on an airplan*^ ^,7^.ng at Hannover (ref- 
erence 27 ) . 

Enri rnd Gruschwitz have now made, in somewhat differ- 
ent manner, the very important attempt to develop purely 
mathematical methods for c-lculnting the co\irse of the r-he- 
nomena in the frictional layer. Buries method iv> simpler, 
while that of Griischwitz is more complete. Lack of space 
forlDidsfurther consider ^'t ion here of these rather compli- 
cated calculations. With these methods it is possi'ble to 
predict the course of the frictional layer for any given 
pressure distrioution and, under some circurastances , even 
to make the important determination as to whether this flow 
will adhere to the wall, as assumed, or will separate at 
some point. A further attempt is now "being made to predict 
in this way the actual characteristics of an airplane wing 
including the profile drag and maximum lift. Should the 
results show a satisfactory agreement with experimental re- 
sults, this method would constitiite a very considera"ble ad- 
vance. 



?URTEEH PROBLEMS 



The investigation of currents in strongly curved 
channels (references 30 and 3l) shows that, aside from the 
"secondary c^irrents" on the side walls as already described 
hy earlier writers , even the" real nature of the turbulence 
is her e substantially altered. The two kinds of phenomena 
are related in that the faster portions of the fluid along 
the curved wall develop stronger centrif\igal forces than 
the slower portions. The faster portions therefore tend to 
displace the slower portions on the outer ;wall. However, 
since the portions in iramediate contact with the wall are 
continiially retarded "by friction, a . materially accelerated 
exchange is produced on the outer, sade of the channel hy 
the displacement of these retarded':port ions . On the con- 
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trary the si owetr' po**r t ions t end toward the inner side and 
the exchange is considerably retarded. 

The phenomena are very similar to those in the flow 
of a .fluid over a heated or cooled bottom surface. In 
the former case the heated and simultaneously retarded 
portions tend to rise from the bottom, while in the latter 
case the' cooled port ions , because of their greater density, 
tend to' remain near the bottom (references 11 and 32), so 
thiat the- turbulent friction is increased in the former 
case and decreased in the latter case. Since both groups 
of phenomena have been or are being investigated in G-ot- 
tingen', numerical -expressions for these influences may be 
exp ec t ed . 

Another important kin>d of phenomena is involved in 
the turbulent spreading of fliiid jets and the wakes of mov- 
ing bodies. The out 'er po rt ions of a jet emerging, e.g. 
from a larger orifice (nozzle, etc.) a-re very unstable and 
develop into a more or less irregular vortex system. Even 
for this kind of phenomena the conception of the mixing 
path held good, and it was possible, with the aid of the 
simple assumption that the mixing path in a cross section 
is constant and proportional to the width of the mixing 
zone at that point,- to predict the form of the mixing zone 
and the velocity distribution in it in a very satisfactory 
manner, whereby only the ratio of the mixing path to the 
mixing 'zone had to be taken from the experiments (refer- 
ences 9, 10, 28, 29, 36) . 

• 

The heat exchange is quite closely related to the 
turbulent velocity exchange. Insofar as it concerns the 
flow along a wall, as shown by the experiments of Elias 
(reference 33), the exchange factor has exactly the same 
value, so that the curve of the temperature distribution 
agrees with the velocity distribution. For the phenomena 
in the wake, of moving bodies, Taylor (r<;fer.ence 34) has 
recently shown that here the heat exchange is twice as 
great as the velocity exchange, so that the temperature 
and velocity curves differ appreciably^. Taylor could also 

^Tf^ylor demonstrates that in this case the. rotational force 
of the main motion is exchanged in the same manner as the 
"heat. The* .exchange factor is p ( du/dy ) ; the rotational 
force in parallel motion, however, is du/dy; and the fall 
of the rotational strength in the direction y is there- 
fore d^^u/dy-^. Taylor shows that then 
(Concluded at bottom of page 23) 
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shoT^ that theoretically the former condition (like form of 
these, curves) is to te expectod when the vortex axes of 
the interference motion are parallel to tne streamlines of 
the main motion, hut the latter (unlii^e) when they are per- 
pendicular to them. The unpublished Crottingen experiments 
of P. Ruden show that the Taylor law of exchan^^e is also 
valid for the spreading of jets. 

It follows therefore that, on closer inspection, there 
are t^'O kinds of turhulence to be distinguished, which dif- 
fer in their nature. We may call one "wall turbulence" and 
the other "jet turbulence." In the former (according to 
Elias) the vortices parallel to the streamlines obviou.sly 
predominate. This rather important discovery will perhaps 
once more indicate the way to a real theory of the phenomena. 
So long as this is not discovered, we must be staisfied with 
half-empirical considerations of the kind here described. 



Translation by D:vi^.^ht M. Miner, 
Nationc.l Advisory Committee 
for Aeronaut ics . 



-i'_-2 = c I'^'i ^lJU which can be integrated to 
3 y " dy dy'^ 

- 0 I" ('—"^for I constant in a cross section. The 

2 Vdy / 

factor ^ in tliis formula differentiates it from our equa- 
tion 2. ^ 
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figure 1. -Velocity profile with turning 
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Pig:ure 2 •-Development of txir- 
"krulence from an in- 
itial distiir'baQce. 



Figure 4.- 
DiBtrilni- 
tion of 
mixing 
path in 
tube at 
large 
Reynolds 
ITumhere. 



Figure 5. -Course "t^jj^ax""*^ ^* large 

Reynolds Hunters. 
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Velocity 

distrilra- 

tion in a 

smooth 

tube. 
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Figure 8.-Resistance in smooth tube according to eq.20 with test values. 




Figure 9. -Resistance coefficients X of rough tuhes. 
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Figure 10. -Roughness function. 



